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ACT I: Space Warps 
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The Challenge of Interstellar Flight 
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VOYAGER 2 


Voyager 1 mission: 

0.7221 spacecraft launched 
in 1977 to study outer solar 
system and boundary with 
interstellar space. 

- After 30+ years, Voyager 1 
is currently at ~120 
Astronomical Units (AU) 
from the sun travelling at 
3.6 AU per year, 

If Voyager 1 were on a I 
trajectory headed to one of 
the Sun's nd3T§?f 
neighboring star systems, 
Alpha Centauri at 4.3 light 
years (or 271,931 AU), it 
would take ~75,000 years to 
traverse this distance at 3.6 
AU/year. 


Pictures courtesy of NASA, 

http://vovager.ipl.nasa.gov/imagesvideo/images/vov 

http://voyager.jpl.nasa.gov/mission/interstellar.html 




DAEDELUS 







Project Daedelus 
sponsored by 
British 

Interplanetary 
Society in 1970's to 
develop robotic 
interstellar probe 
capable of reaching 
Barnard's star, at ~6 
light years away, in 
50 years. 

The resulting 
spacecraft was 
54,000t, 

92% fuel for fusion 
propulsion system. 

ISS is ~450t 



Picture courtesy of Adrian Mann, www.bisbos.com. used with permission 
ISS picture courtesy of NASA (STS-135 fly-around) 




% - * ♦ 

T ~T Tf 7T 












* 



rx 

NASA 



Hyper-fast interstellar travel 



• Is there a way within the framework of physics such that one 
could cross any given cosmic distance in an arbitrarily short 
period of time, while never locally breaking the speed of light 
(11 th commandment)? 

SPACEWARPS 
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Images courtesy NASA 
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Inflation: Alcubierre Metric 1 
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Warp Metric: 




£9 


3 * If 

«cO’' y> 


Appealing Characteristics 


Proper acceleration 
in the bubble is 
formally zero 


Images courtesy NASA 



Flat space-time inside the bubble 


Unappealing 
characteristic 

(square peg, round hole) 

MCC clocks synchronized with onboard 

clocks 


(divergence of phi = 0) 


(Coordinate time = proper time) 




"bubble" thickness decreases 



"bubble" thickness decreases 

10 

Surface plots of York Time & T 00 , <v>=10c, 10 meter diameter volume, variable warp "bubble" thickness 
















































As bubble thickness increases, 
York Time intensity decreases 



Allowing the bubble to get 
thicker reduces the flat space- 
time real-estate in the center 


Changing topology greatly 
reduces the energy required 



But space-time is really stiff: c*/87tG 

Can we further reduce the energy 
required by reducing the stiffness? 


Maybe... but we need to engage 
higher dimensional models to do so 
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Bubble Oscillation Optimization 
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Oscillate the bubble intensity 
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Exotic Mass Warp Requirements, 10m diameter, v =10c 
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Topology 


Approximations 
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Jupiter 1.9X10 37 
Earth 6.0x10“ 

Moon 7.4X10 22 


Bulk Velocity, dU/dt 

■ dU/dt =0c 

* dU/dt =0.9c 

x dU/dt =0.99 9 c 
-*■ dU/dt =0.99 999c 
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Nimitz carrier 9.lxlO 7 


ISS 4.0xl0 5 


a 


Voyager 



** A 


| 

L 

/ ' 

13\ 











































































































Cv 

NASA 




Warp Field Interferometer 




mr 


He-Ne Laser 
~ 633 nm 


toroid capacitor ring 
normal to beam path 

resultant spherical field 



beam splitter 


detector 





Warp Field Interferometer developed after 
putting metric into canonical form 1 : 


ds 2 =(v s 2 f(r s ) 2 -\) 


dt- 


v s f(r s ) 


v;/(0 2 -1 


dx 


2 2 7 

— dx + dy + dz 


Generate microscopic warp bubble that 
perturbs optical index by 1 part in 
10,000,000 

Induce relative phase shift between split 
beams that should be detectable. 



Numerical simulation 
of fringe at detector 
with device off 


Numerical simulation 
of fringe at detector 
with device on 


2D Analytic 
Signal 
processing 
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M(x, y)e 


Numerical simulation 
of contour plot of (j) 
for warp field 


White, H., "A Discussion on space-time metric engineering/' Gen. Rel. Grav. 35, 2025-2033 (2003). 















2 Polarizers 

(jointly used for 
intensity 
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Fabry-Perot 


Consists of two reflecting, 
highly parallel surfaces, called 
an Etalon 

The interference pattern is 
created within the Etalon 

Multiple reflections in the 
Etalon reinforce the areas 
where constructive and 
destructive interference 
occurs 

Allows for much higher- 
precision measurements of 
fringes (image averaging 
without software) 


Example: Michelson-Morley 
Interferometer image for 
Sodium source 


Example: Fabry-Perot 
Interferometer image for 
Sodium source (note doublet) 
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Frequency Spectrum (363,369) 


Frequency (Hz) 




















Primary/Auxiliary 


Isolator pressure 


Height control pressure 
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FFT of imager data at frequency of i 













n*M Open-air etalon Implementation 





Frequencies 
of interest 




Inactive 


Active 
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Time of Flight Schematic 



200MHz test device 



Detectors 


Thorlabs High-Speed Avalanche Detector 



Model 

APD210 

Rise Time 

0.5 ns 

Supply Voltage 

+12 to +15 V 

Current Consumption 

200 mA 

Max. Incident Power 

10 mW 

Spectral Range 

400-1000 nm 

Frequency Range 

1-1600 mHz 

Maximum Gain 

2.5xl0 5 V/W 


HeNe 


EO Modulator 

Newport Broadband Amplitude Modulator 


Model 

Type 

Operating Frequency 


4102NF 

Broadband Amplitude Modulator 
DC-200 MHz 


Wavelength Range 
Material 
Maximum Vn 


500-900 nm 
MgO:LiNbO s 
195 V @ 633 nm 


Maximum Input Power 2 W/mm 2 @ 532 nm 


Aperture Diameter 

2 mm 

RF Bandwidth 

200 MHz 

RF Connector 

SMA 

Input Impedance 

10 pF 

Maximum RF Power 

10W 

Connector 

SMA 




Agilent Technologies Infiniium DS09254A 2.5 GHz Oscilloscope 

•2.5 GHz bandwidth across all 4 analog channels 
•20 GSa/s max. sample rate 

•Standard 20 Mpts memory per channel, upgradeable to 1 Gpts 
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Forward Plan 


NASA 


Explore the d<fi/dt dependency in future test 
devices 

- The idea of an optimized space warp needs vacuum 
energy, and large d(j)/dt 

— both of these conditions are predicted to be present in 
the q-thruster technology also being explored in the 
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Original Matthew Jeffries concept from mid 
1960's, rendered by Mark Rademaker 
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RademakerXdrex files 


Matthew Jeffries is the artist that created the familiar Star Trek enterprise 
look 
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on Dr. White's theoretical findings, rendered by 
with artwork and inputs from Mike Okuda 
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Act II: 

Experimental Thrust Measurement 

of RF Test Articles 
exploring q-thruster models 



Q-Thruster Background 

> A Q-thruster is a form of electric propulsion 

> Through the use of electric and magnetic fields, a Q-thruster pushes quantum 
particles (electrons/positrons) in one direction, while the Q-thruster recoils to 
conserve momentum 

> This principle is similar to how a submarine uses its propeller to push water 
in one direction, while the submarine recoils to conserve momentum 

> Based on test and theoretical model development, expected thrust to 
power for initial flight applications is 0.4N/kW (~7x Hall) 

> 0.4 N/kW enables power-constrained HEO SEP missions to close without 
needing chemical kick stages and very long transit times. 

> 0.4N/kW coupled with persistent power (e.g. NEP) enables rapid transit 
missions throughout the solar system. 


"If vacuum fluctuations can be harnessed for propulsion by anyone besides 
science-fiction writers, the purely engineering problems of interstellar flight 
would be solved." - A.C. Clarke 






Magnetic Damper 


Electrostatic Fins 





down to 5xlO A -6 Torr 


Vacuum 


Rackmount DAQ 


Linear Displacement Sensor & Adjustment Remote 

























Cannae RF Test Article 

> Test article is a pillbox/beam pipe design 
fashioned after an RF resonant cavity design 
used in high energy particle accelerators. 

> Each test article is -11" in diameter and 4-5" 
between ends of the beam pipes, not counting 
beam pipe extensions or antenna mounts. 

> Concept is a modified particle accelerator 
design that incorporates engraved radial slots 
along the outside edge of the resonant cavity 
interior, but on only one side of the pillbox 
(equatorially-asymmetric). 

> Cannae theorizes that asymmetric engraved 
slots would result in a force imbalance 
(thrust). 

> As a result, a second (control) test article was 
fabricated without the internal slotting 
dubbed the null test article. 
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Testing 

Overview 



> RF signal fed into one of test 
article beam pipes via an 
approximately 20-gauge copper 
wire drive antenna. 

> On opposite beam pipe, an 
approximately 20-gauge copper 
wire serves as a sense antenna 
for manual tuning during 
testing. 


> The longer beam pipe is the RF drive 
antenna that in practice ends up being 
a l A wave resonance system in its own 
right and has a dielectric PTFE slug in 
the throat in both the slotted and null 
test article. 

> It is this characteristic that became an 
item of further consideration after 
completion of the test campaign. 




ZX73-2500M-S+ 
10-to-2,500 MHz 
W-Atten. 


ZHL-100W-13+ Temp-1 
800-to-1,000 MHz 
100W (50dBm) Amp / 


Spectrum 

Analyzer 


Cannae Fwd. / Rev. 
Directional 


RF out 


ZX47-40LN S+ 
10-8,000 MHz 
Power Detector # 1 

/b 

ZX47 40LN-S+ 
10-8,000 MHz 
Power Detector #2 


(Stainless Steel Chamber) 


Given: 1U, LM Contacts & #22 
AWG wiring implies 7.0 Amps 
max per chassis circuit. 
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Representative Test Run 


















































































r 85.6 


73.6 
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-85 0 
-845 
-840 
-83.5 
-83.0 
-82.5 
-82.0 
-81 5 
-81.0 
-80.5 
-80 0 
-79.5 
-79.0 
-78.5 
-78.0 
-77.5 
-77.0 
-76.5 
-76 0 
-75.5 
-75.0 
-74.5 
-74.0 
-73.5 
-73.0 
-72.5 


L 71.9 



Test article does not get 
warm during operation 
(high Q resonance system) 
Amplifier kept below 100 F 







Cannae Test Results Summary 

> The resistive RF Load evaluation indicated no significant 
systemic cause for torsion pendulum displacement. 

> Based upon this observation, both test articles (slotted and 
unslotted) produced significant thrust in both orientations 
(forward and reverse). 

> Test schedule constraints prevented multiple data points 
to be gathered in the reverse orientation, and the single 
data point for each test article is insufficient to allow 
comparative conclusions (between slotted and unslotted) 
to be drawn. 

> However, for the forward thrust orientation, the difference 
in mean thrust between the slotted and unslotted was less 
than two percent. 

> Thrust production did not appear to be dependent on the 
slotting. 



Configuration 1A 



Configuration 

Test Article 

Thrust 

Direction 

Thrust Range 
(f^N) 

Mean 

Thrust 

(uN) 

Number of 

Test Runs 

1A 

Slotted 

Forward 

31.7-45.3 

40.0 

5 

IB 

Slotted 

Reverse 

48.5 

48.5 

1 

2A 

Unslotted 

Forward 

35.3 - 50.1 

40.7 

4 

2B 

Unslotted 

Reverse 

22.5 

22.5 

1 

RF Load 

50Q Load 

N/A 

0.00 

0.00 

2 


Configuration 2B 



RF Dummy Load 
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COMSOL Analysis of Cannae Test Articles 


> 


Computer modeling of the electric field 
within the pillbox and beam pipe 
illustrates weakness of electric field in 
vicinity of cavity slots and relative strength 
of the electric field within the beam pipe, 
especially in the drive antenna coaxial 
cable and the region around the cable 
within the PFTE dielectric slug as seen in 
the following figure. 



> Consideration of the dynamic fields in the wave resonance tube shows that there 
is always a net Poynting vector meaning that the RF launcher tube assembly with 
dielectric cylinder common to both the slotted and smooth test articles was 
potentially a Q-thruster where the pillbox is simply a matching network. 


> The predicted thrust using the Q-thruster analytic model is 34 micronewtons with 
input power of 25 Watts and Quality factor of 8000. 
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Tapered RF Q-thruster 

> Just like there are many ways to build a race car motor, there are 
many ways to build a thruster, or in this case, a Q-thruster. 

> The figures depict one of the early COMSOL models representing 
an early possible construction of a tapered RF unit alongside the 
actual construction that was finally implemented as informed by 
the COMSOL findings. 

> The tapered RF geometry was explored as it pertains to some recent 
findings published by the Northwest Poly technical University. 
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4" dielectric 
resonator 



TE012 mode Test article 
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Systemic Forces on the Pendulum-4 
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Distance and Chart Data 
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581.000 
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Thermal response 
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29.1uN 


Browse 
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576.000- 
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200V, 

29.1uN 








File Name 


Points Saved 0 
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COMSOL RF Analysis 














TM212,1,937.115MHz, ~30W, Q L =6,726, 


Microns Distance and Chart Data 































































































































































































TM212,1,937.115MHz, ~35W, Q L =6,726, 
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Forward Thrust Campaign 






















FORWARD WORK 

> JSC: Complete the vacuum test campaign 


> Perform an IV&V test campaign at the Glenn Research Center using 
their low thrust torsion pendulum followed by a possible repeat 
campaign at the Jet Propulsion Laboratory using their low thrust 
torsion pendulum. 


> The Johns Hopkins University Applied Physics Laboratory has also 
expressed an interest in performing a Cavendish Balance style test 
with the IV&V shipset. 
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GRC Torsion Pendulum 






Act III: 

Value Proposition 



■ Dawn propulsion & 
power system 

■ Three NSTAR xenon 
ion thrusters 

■ 0.04 N/kWe thrust- 
to-power ratio 

■ ~90 mN/thruster 

■ 2.3 kWe solar power 
at 1 AU 

■ 425 kg of xenon 
propellant 

Dawn with 2.3kWe Ion Thrusters Dawn with 2.3kWe O-thrusters 



Activity 

Dates 

Segment 

(months 

) 

Elapsed 

(months) 

Dates 

Segment 

(months) 

Elapsed 

(months) 

Earth to 
Vesta 

Sep 2007- 
Jul 2011 

46 

46 

Sep 2008 - 
Apr 2009 

7 

(15% of ion) 

7 

Vesta orbit 

Jul 2011 - 
Jul 2012 

12 

(fixed) 

58 

Apr 2009 - 
Apr 2010 

12 

(fixed) 

19 

Vesta to 
Ceres 

Jul 2012- 
Feb 2015 

31 

89 

Apr 2010 - 
Jul 2011 

15 

(48% of ion) 

34 

Ceres orbit 

Feb-Jul 2015 

5 

(fixed) 

94 

(7.8 yrs) 

Jul-Dec 2011 

5 

(fixed) 

39 

(3.3 yrs) (41%) 
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300 kW SEP Mars 


• 70t stack departs from 
DRO 

• 300kW SEP 

• 0.4N/kW Q-thrusters 

• 50-day stay in Deimos 
orbit around Mars 

• Total mission duration of 
788 days and 2AU 
maximum distance from 
the sun. 
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2MW NEP Mars 


JrtARS 

66 days 


ARTH 



90t spacecraft 
2MW power 
0.4N/kW (800N) 

246 day mission with 
70 day stay at Mars 



90t spacecraft: 

• 50t cargo 

• 20t 2MW nuclear reactor (10 kg/kW) 

• 20t Q-thruster bank (lOkg/kW) 
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Crewed Titan/Enceladus Mission 

• Same vehicle characteristics as Jupiter mission assumed 

• Trip time to/from Saturn ~9 months each way, 6 month stay in Europa 
vicinity, 6 month stay in Enceladus vicinity. 32-month total mission. 


Assuming: 

T s =0.4 N/kWe 
P=2000 kWe 

°Wr = 10 kg/kWe 
a prop = 10 kg/kWe 

m pi = 50t 


Shorter than "current" conjunction-class Mars missions 



T/m =0.91 m-g's 
Vehicle mass = 901 











Act IV: 

Dynamics of the Vacuum and Casimir 
Analogs to the Hydrogen Atom 



Principles of Q-thrustei 


> Principle 1 : Local mass concentrations, say in the 
form of a conventional capacitor with a ceramic 
dielectric, affect vacuum fluctuation density 
according to equation 1 


> Principle 2 : Just as relativistic acceleration (Unruh 
radiation) can change the apparent relative density 
of the vacuum, so too can higher order derivatives 
according to equation 2. 


> Principle 3 : The tools of MagnetoHydroDynamics 
(MHD) can be used to model this modified 
vacuum fluctuation density analogous to how 
conventional forms of electric propulsion model 
propellant behavior. 
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TABLE I: This table shows the derived “density” of a 


given energy state n, with Z — 1 


n a 

radms(m) 

E(eV) 


E(J) 


P {kd/ 

m 3 ; 

) 

1 

5. 

29 

X 

10“ 

TI 

13.60 

2. 

.176 

X 

10“ 

-18 

3.905 

X 

10' 

-5 

2 

2. 

11 

X 

10“ 

■10 

3.40 

5. 

.440 

X 

10“ 

-19 

1.525 

X 

10' 

-7 

3 

4. 

76 

X 

10“ 

-10 

1.51 

2. 

.418 

X 

10“ 

-19 

5.952 

X 

10- 

-9 

4 

8. 

46 

X 

10“ 

■10 

0.85 

1. 

.360 

X 

10“ 

-19 

5.959 

X 

10“ 

10 

5 

1 

.32 

X 

10' 

-9 

0.54 

8. 

.704 

X 

10“ 

-20 

9.997 

X 

10“ 

11 

6 

1 

.90 

X 

10‘ 

-9 

0.38 

6. 

.044 

X 

10“ 

-20 

2.325 

X 

10“ 

11 

7 

2 

.59 

X 

10' 

-9 

0.28 

4. 

.441 

X 

10“ 

-20 

6.774 

X 

10“ 

■12 


a The primary quantum number n is only varied from 1 to 7 here. 


47reo n 2 h 2 
e 2 m e 

= n 2 5.29 x 1 CP 11 meters , n = 1,2,3,... 


E n — — 


1 


m 

( e2 VI 

1 

2h 

\47re 0 y 

n 2 


=-^13.6eV, n— 1,2,3, 




1 Z.n 



5.29x10 11 meters , n = 1,2,3,... 



Z 2 

—s-13.6eV, n = 1. 2,3,... 
n 1 




n 


9 4 

C 3^ 


.3 
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TABLE II: This table compares the derived “density” 
of a given energy state n, with Z = 1 to the Casimir 
density for a cavity with a separation distance of 2r n . 


n 

radius (m) 

p (kg/m 3 

) Casimir {kgj 

yr 

Eatio h 

1 

5.29 

X 

10“ 

7 iL 

3.91 

X 

10 - 

-5 

1.16 

X 

10 - 

-4 

2.96 

2 

2.11 

X 

10“ 

10 

1.53 

X 

10 - 

-7 

4.51 

X 

10 - 

-7 

2.96 

3 

4.76 

X 

10“ 

10 

5.95 

X 

10 - 

-9 

1.76 

X 

10 - 

-8 

2.96 

4 

8.46 

X 

10“ 

10 

5.96 

X 

10“ 

10 

1.76 

X 

10 - 

-9 

2.96 

5 

1.32 

X 

10 ' 

-9 

1.00 

X 

10“ 

10 

2.96 

X 

10 “ 

10 

2.96 

6 

1.90 

X 

10' 

-9 

2.33 

X 

10“ 

■11 

6.88 

X 

10 “ 

11 

2.96 

7 

2.59 

X 

10 ' 

-9 

6.77 

X 

10“ 

12 

2.00 

X 

10 “ 

11 

2.96 


a This is the Casimir force per unit area multiplied by 1/cr . 

This ratio is the Casimir column value divided by the p column 
value 


Friedmann equation: 


a 

a 


4t tG 


3 


(pc 2 + 3 P) 


P < — /9C 2 /3 
w ~ —1/3 
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Z,n 


1 hen 2 


3c 2 240d 4 





































































Vacuum density perturbation (kg/m A 3) 
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Mesh sizes: 

1.5 



n=l: 1 pm 
n=2: 2.5 pm 
n=3: 5 pm 

ij 



n=4: 10 pm 
n=5: 25 pm 

0.5" 


t \ •' i 

n=6: 50 pm 
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-0.5" 
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FIG. 6: COMSOL 2D axisymmetric model: element 
sizes for n — 1 to n = 7 is 1 pro, 2.5pm, 5pro, 1 Opro, 
25pro, 50pro, and 50pm respectively. In this figure, the 
mesh size is too dense to be discernable for n = 5 and 
lower. FIG. 4 shows the mesh for the n — 1 and n = 2 


regions 



FIG. 4: Close-up of COMSOL 2D axisymmetric model. 


Eigenfrequency = S.394401el4 




field (Pa) 

▲ 3 . 981 . 


3 . 981 . 


KHI 


3 . 423 ! 


3 145 ! 




torn 


(a) Macroview 


(b) Microview 


FIG. 5: COMSOL analysis results for n = 1 
eigenfrequency: panel 5a shows the model out to the 
n = 6 orbital, and panel 5b show's a close-up view of the 

n = 1 solution 

























































TABLE IV: This table shows the predicted “orbital” 
frequency and the COMSOL analysis eigenfrequend.es 
for the n = 1 to n = 7 states for hydrogen. 


n 

Orbital /reg a 

COMSOL freq h 

%error 

i 

6.58 

X 

10 15 

6.25 

X 

10 1 * 

-4.98 

2 

8.23 

X 

10 14 

8.23 

X 

10 14 

0.05 

3 

2.44 

X 

10 14 

2.38 

X 

10 14 

-2.48 

4 

1.03 

X 

10 14 

1.01 

X 

10 14 

-1.59 

5 

5.26 

X 

10 13 

4.98 

X 

10 13 

-5.36 

6 

3.05 

X 

10 13 

3.48 

X 

10 13 

14.28 

7 

1.92 

X 

10 13 

2.13 

X 

10 13 

11.16 


a Orbital frequency is in Hz. 
b COMSOL frequency is in Hz. 
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Eigenfreguency -A 982496*13 Eigenfreguency-3.481617*13 Dgenfr*quency-2.132714el3 91755.6094Bli 

jrface; Total acoustic pressure field (Pal Surface: Total acoustic pressure field (Pal Surface: Total acoustic pressure field (Pa) 



FIG. 7: COMSOL analysis results of the acoustic 
“natural” vacuum model The orbital shells (dark lines) 
can be counted, but the n = 1 radius is quite small as 
seen from the top left thumbnail that depicts the 
COMSOL eigenfrequency solution for that orbital. 
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FIG. 9: 2D Axisymmetric model results that capture 
axisymmetric acoustic solutions like 2 p, 3d, and 4/ 

orbitals (m = 0). 



FIG. 3: Plot for the Z = 1, 2p orbital from Orbital 
Viewer Software. 

































